function Zhu2022GAD_combined_Cell_Recurring_Prob(Data_Summary)
%% get all basic info to put into a summary table ready
animalNames = {'G2_F4_18'; 'G2_F5(2)_28'; 'G2_F5(2)_31'; 'G2_F5(2)_34'; 'G2_F6(2)_42'; 'CCK_F6_37'; 'CaMKIIa_F6_32'; 'G2_F4_17'; 'G2_F5(2)_21'; 'G2_F5(2)_45'; 'G2_F5(2)_53'; 'G2_F6(2)_36'; 'CaMKIIa_F6_19'; 'CaMKIIa_F6_20'; 'CaMKIIa_F6_21'; 'CaMKIIa_F6_22';};
genotype = {'p-WT'; 'p-WT'; 'p-WT'; 'p-WT'; 'p-WT'; 'p-WT'; 'p-WT'; 'a5-i-KO'; 'a5-i-KO'; 'a5-i-KO'; 'a5-i-KO'; 'a5-i-KO'; 'a5-pyr-KO'; 'a5-pyr-KO'; 'a5-pyr-KO'; 'a5-pyr-KO'};

%% collect saline data
[G18_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'saline');
[G28_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'saline');
[G31_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'saline');
[G34_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'saline');
[G42_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'saline');
[CCK37_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'saline');
[CaMKIIa32_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'saline');
[G17_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'saline');
[G21_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'saline');
[G45_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'saline');
[G53_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'saline');
[G36_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'saline');
[CaMKIIa19_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'saline');
[CaMKIIa20_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'saline');
[CaMKIIa21_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'saline');
[CaMKIIa22_saline_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'saline');

% concatenate results
saline = cat(1, G18_saline_M, G28_saline_M, G31_saline_M, G34_saline_M, ...
    G42_saline_M, CCK37_saline_M, CaMKIIa32_saline_M, G17_saline_M, G21_saline_M, ...
    G45_saline_M, G53_saline_M, G36_saline_M, CaMKIIa19_saline_M, CaMKIIa20_saline_M, ...
    CaMKIIa21_saline_M, CaMKIIa22_saline_M);

%% collect etomidate 2mg/kg data
[G18_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'etomidate 2mg/kg');
[G28_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'etomidate 2mg/kg');
[G31_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'etomidate 2mg/kg');
[G34_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'etomidate 2mg/kg');
[G42_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'etomidate 2mg/kg');
[CCK37_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'etomidate 2mg/kg');
[CaMKIIa32_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'etomidate 2mg/kg');
[G17_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'etomidate 2mg/kg');
[G21_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'etomidate 2mg/kg');
[G45_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'etomidate 2mg/kg');
[G53_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'etomidate 2mg/kg');
[G36_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'etomidate 2mg/kg');
[CaMKIIa19_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'etomidate 2mg/kg');
[CaMKIIa20_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'etomidate 2mg/kg');
[CaMKIIa21_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'etomidate 2mg/kg');
[CaMKIIa22_etom2_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'etomidate 2mg/kg');

% concatenate results
etom2 = cat(1, G18_etom2_M, G28_etom2_M, G31_etom2_M, G34_etom2_M, ...
    G42_etom2_M, CCK37_etom2_M, CaMKIIa32_etom2_M, G17_etom2_M, G21_etom2_M, ...
    G45_etom2_M, G53_etom2_M, G36_etom2_M, CaMKIIa19_etom2_M, CaMKIIa20_etom2_M, ...
    CaMKIIa21_etom2_M, CaMKIIa22_etom2_M);

%% collect etomidate 4mg/kg data
[G18_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'etomidate 4mg/kg');
[G28_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'etomidate 4mg/kg');
[G31_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'etomidate 4mg/kg');
[G34_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'etomidate 4mg/kg');
[G42_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'etomidate 4mg/kg');
[CCK37_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'etomidate 4mg/kg');
[CaMKIIa32_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'etomidate 4mg/kg');
[G17_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'etomidate 4mg/kg');
[G21_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'etomidate 4mg/kg');
[G45_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'etomidate 4mg/kg');
[G53_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'etomidate 4mg/kg');
[G36_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'etomidate 4mg/kg');
[CaMKIIa19_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'etomidate 4mg/kg');
[CaMKIIa20_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'etomidate 4mg/kg');
[CaMKIIa21_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'etomidate 4mg/kg');
[CaMKIIa22_etom4_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'etomidate 4mg/kg');

% concatenate results
etom4 = cat(1, G18_etom4_M, G28_etom4_M, G31_etom4_M, G34_etom4_M, ...
    G42_etom4_M, CCK37_etom4_M, CaMKIIa32_etom4_M, G17_etom4_M, G21_etom4_M, ...
    G45_etom4_M, G53_etom4_M, G36_etom4_M, CaMKIIa19_etom4_M, CaMKIIa20_etom4_M, ...
    CaMKIIa21_etom4_M, CaMKIIa22_etom4_M);

%% collect etomidate 6mg/kg data
[G18_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'etomidate 6mg/kg');
[G28_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'etomidate 6mg/kg');
[G31_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'etomidate 6mg/kg');
[G34_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'etomidate 6mg/kg');
[G42_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'etomidate 6mg/kg');
[CCK37_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'etomidate 6mg/kg');
[CaMKIIa32_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'etomidate 6mg/kg');
[G17_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'etomidate 6mg/kg');
[G21_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'etomidate 6mg/kg');
[G45_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'etomidate 6mg/kg');
[G53_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'etomidate 6mg/kg');
[G36_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'etomidate 6mg/kg');
[CaMKIIa19_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'etomidate 6mg/kg');
[CaMKIIa20_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'etomidate 6mg/kg');
[CaMKIIa21_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'etomidate 6mg/kg');
[CaMKIIa22_etom6_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'etomidate 6mg/kg');

% concatenate results
etom6 = cat(1, G18_etom6_M, G28_etom6_M, G31_etom6_M, G34_etom6_M, ...
    G42_etom6_M, CCK37_etom6_M, CaMKIIa32_etom6_M, G17_etom6_M, G21_etom6_M, ...
    G45_etom6_M, G53_etom6_M, G36_etom6_M, CaMKIIa19_etom6_M, CaMKIIa20_etom6_M, ...
    CaMKIIa21_etom6_M, CaMKIIa22_etom6_M);

%% collect etomidate 7mg/kg data
[G18_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'etomidate 7mg/kg');
[G28_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'etomidate 7mg/kg');
[G31_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'etomidate 7mg/kg');
[G34_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'etomidate 7mg/kg');
[G42_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'etomidate 7mg/kg');
[CCK37_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'etomidate 7mg/kg');
[CaMKIIa32_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'etomidate 7mg/kg');
[G17_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'etomidate 7mg/kg');
[G21_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'etomidate 7mg/kg');
[G45_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'etomidate 7mg/kg');
[G53_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'etomidate 7mg/kg');
[G36_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'etomidate 7mg/kg');
[CaMKIIa19_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'etomidate 7mg/kg');
[CaMKIIa20_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'etomidate 7mg/kg');
[CaMKIIa21_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'etomidate 7mg/kg');
[CaMKIIa22_etom7_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'etomidate 7mg/kg');

% concatenate results
etom7 = cat(1, G18_etom7_M, G28_etom7_M, G31_etom7_M, G34_etom7_M, ...
    G42_etom7_M, CCK37_etom7_M, CaMKIIa32_etom7_M, G17_etom7_M, G21_etom7_M, ...
    G45_etom7_M, G53_etom7_M, G36_etom7_M, CaMKIIa19_etom7_M, CaMKIIa20_etom7_M, ...
    CaMKIIa21_etom7_M, CaMKIIa22_etom7_M);

%% collect etomidate 8mg/kg data
[G18_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'etomidate 8mg/kg');
[G28_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'etomidate 8mg/kg');
[G31_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'etomidate 8mg/kg');
[G34_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'etomidate 8mg/kg');
[G42_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'etomidate 8mg/kg');
[CCK37_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'etomidate 8mg/kg');
[CaMKIIa32_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'etomidate 8mg/kg');
[G17_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'etomidate 8mg/kg');
[G21_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'etomidate 8mg/kg');
[G45_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'etomidate 8mg/kg');
[G53_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'etomidate 8mg/kg');
[G36_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'etomidate 8mg/kg');
[CaMKIIa19_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'etomidate 8mg/kg');
[CaMKIIa20_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'etomidate 8mg/kg');
[CaMKIIa21_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'etomidate 8mg/kg');
[CaMKIIa22_etom8_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'etomidate 8mg/kg');

% concatenate results
etom8 = cat(1, G18_etom8_M, G28_etom8_M, G31_etom8_M, G34_etom8_M, ...
    G42_etom8_M, CCK37_etom8_M, CaMKIIa32_etom8_M, G17_etom8_M, G21_etom8_M, ...
    G45_etom8_M, G53_etom8_M, G36_etom8_M, CaMKIIa19_etom8_M, CaMKIIa20_etom8_M, ...
    CaMKIIa21_etom8_M, CaMKIIa22_etom8_M);

%% collect altered contexts data
[G18_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_18', 'altered contexts');
[G28_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_28', 'altered contexts');
[G31_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_31', 'altered contexts');
[G34_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_34', 'altered contexts');
[G42_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_42', 'altered contexts');
[CCK37_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CCK_F6_37', 'altered contexts');
[CaMKIIa32_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_32', 'altered contexts');
[G17_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F4_17', 'altered contexts');
[G21_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_21', 'altered contexts');
[G45_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_45', 'altered contexts');
[G53_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F5(2)_53', 'altered contexts');
[G36_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'G2_F6(2)_36', 'altered contexts');
[CaMKIIa19_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_19', 'altered contexts');
[CaMKIIa20_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_20', 'altered contexts');
[CaMKIIa21_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_21', 'altered contexts');
[CaMKIIa22_altered_M] = Cell_Recur_Prob_Mean_Prism(Data_Summary, 'CaMKIIa_F6_22', 'altered contexts');

% concatenate results
altered_contexts = cat(1, G18_altered_M, G28_altered_M, G31_altered_M, G34_altered_M, ...
    G42_altered_M, CCK37_altered_M, CaMKIIa32_altered_M, G17_altered_M, G21_altered_M, ...
    G45_altered_M, G53_altered_M, G36_altered_M, CaMKIIa19_altered_M, CaMKIIa20_altered_M, ...
    CaMKIIa21_altered_M, CaMKIIa22_altered_M);

%% Now, make summary table

Cell_Recur_Prob_Summary_Table = table(animalNames, genotype, saline, etom2, etom4, etom6, etom7, etom8, altered_contexts);

%% Export table

[ds, ~] = getdsap;
path = fullfile(ds.metadataPath, 'Prism excel data', 'Cell Recurring Probability.xlsx');

if ~exist((fullfile(ds.metadataPath, 'Prism excel data')), 'dir')
    mkdir(fullfile(ds.metadataPath, 'Prism excel data'));
end
writetable(Cell_Recur_Prob_Summary_Table, path);

end